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ABSTRACT 

Phenetic affinities of Australasian oystercatchers were eluci¬ 
dated by multivariate statistical analysis of seven morpho¬ 
metric characters taken from museum skins. Nineteen opera¬ 
tional taxonomic units (OTUs), representing all the currently 
recognized species and most of the subordinate taxa, were 
used to calibrate morphological variation in Australasian 
Haematopus against the range displayed by the Haematopo- 
didae. The rather homogeneous nature of variation among 
OTUs was not well suited for analysis by hierarchical 
clustering methods, but a non-hierarchical ordination method 
utilizing both principal components and nonmetric scaling 
produced excellent summaries of the similarity matrices. 

The South Island Pied Oystercatcher (H. o. finschi ) of 
New Zealand clustered tightly with the Eurasian H. ostra- 
legus and thus seems to be of Palaearctic origin. All the 
remaining Australasian taxa grouped in another cluster with 
the New World forms, and it appears likely that the 
Australian Pied Oystercatcher will have to be split from the 
ostralegus group as a separate species. Comments on the 
taxonomy and zoogeography of the family are postponed 
until a current study of a much larger data set is completed. 

“It cannot be but mischievous to any study to have divers 
systems of nomenclature simultaneously co-existent in the 
same dominion.” Milligan (1911). 

INTRODUCTION 

Australasian oystercatchers are widely distributed 
around the sea shores of New Zealand, Australia, 
most islands in Torres Straits, Misima, Kai and 
Aru Islands, and the southern coast of New Guinea 
(Condon, 1975). The affinities and hence the 
number of taxa that should be recognized from 
this region have been subject to lively debate in 
the past, and remain largely unresolved today. 

I have outlined previously the chaotic nomen- 
datural history of the New Zealand taxa, and have 
suggested that this instability was due primarily 
to an essentially typological assessment of varia¬ 


tion in morphological characters (Baker, 1975). 
The Australian taxa have suffered similar treat¬ 
ment, as illustrated below. The first oystercatcher 
to receive formal recognition was a pied bird, 
described by Vieillot (1817) a s Haematopus Iongi- 
rostris. Three further pied birds from northwestern 
Australia, New South Wales and Queensland were 
described respectively as H. picatus King, 1826, 

H. australianus Gould, 1838 and H. longirostris 
mattingleyi Mathews, 1912. Mathews (1912) syn- 
onymized H. australianus with the nominate H. 

I. longirostris, and restricted picatus as a sub¬ 
species of H. longirostris. The following year he 
submerged these taxa under H. ostralegus longi¬ 
rostris Vieillot, ostensibly because of the paucity 
of material then available by which good sub¬ 
species could be recognized (Mathews, 1913). All 
Australian pied birds were subsequently lumped 
under H. ostralegus Linnaeus, 1758, by Mathews 
& Iredale (1921), and no subordinate taxa were 
designated. Mathews (1927) assigned longirostris 
and picatus as subspecies of H. ostralegus, but 
Peters (1934) considered only H. o. longirostris 
as valid. Oliver (1955) concluded that specimens 
from New South Wales, northern New Zealand 
and the Chatham Islands were all referable to one 
species, H. longirostris. Condon (1975) relegated 
longirostris to a subspecies of H. ostralegus, but 
noted that “the taxonomy is very confused and 
some authors are inclined to follow Oliver”. 

Similar taxonomic confusion is apparent for the 
melanistic oystercatchers of Australia. The original 
description of the Sooty Oystercatcher was given 
by Gould (1845) as Haematopus fuliginosus. 
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Castelnau & Ramsay (1877) erected a new species, 
H. ophthalmicus (lapsus calumni), for a melanistic 
specimen with a large fleshy patch around the eye 
from Bountiful Island in the Gulf of Carpentaria. 
In 1912 Mathews referred both fuliginosus and 
ophthalmicus to subspecies of the New Zealand 
oystercatcher H. unicolor Forster, 1844 and erected 
a new subspecies H. u. bernieri for a specimen 
from Bernier Island in Western Australia. On the 
basis of supposed priority he synonymized unicolor 
with niger in 1913, but retained the above sub¬ 
species. Mathews & Iredale (1921) reversed this 
decision and submerged all taxa under H. uni¬ 
color. In 1927 Mathews reinstated fuliginosus, 
ophthalmicus and bernieri as subspecies of H. 
unicolor, but Peters (1934) rejected unicolor as a 
species epithet and recognized only H. f. fuligi¬ 
nosus and H. f. ophthalmicus from Australia. 
These latter taxa have endured since and are used 
by Condon (1975). 

The cluttering of the literature with numerous 
nomenclatural changes is perhaps understandable 
in the context of the nature and relative recency 
of ornithological discovery in Australasia. How¬ 
ever, if progress is to be made beyond the point 
where taxonomic decisions reflect little more than 
the personal preferences of their proponents, it is 
clear that a substantive study of the systematic 
relationships of the taxa will have to be made. 
The taxonomic instability outlined above stems 
principally from an inadequate assessment of the 
affinities of allopatric taxa, especially those from 
the Old and New Worlds. This paper presents 
a multivariate analysis of the phenetic affinities of 
Australasian oystercatchers, based on external mor¬ 
phometric skin characters traditionally used by 
avian systematists. 

MATERIALS AND METHODS 

Selection of taxa for comparative study 

Because the phenetic affinities of Australasian 
Haematopus need to be cast in the broader per¬ 
spective of the range of morphological variation 
in the Haematopodidae, the selection of taxa for 
comparative study is most important. Where mul¬ 
tivariate statistical methods are used, this sampling 
factor becomes critical because the configuration 
of the multidimensional space defined by the 


characters and/or the taxa is dependent on the 
range of variation spanned. If the placement of 
taxa in such a space is to be a reliable and repeat- 
able representation of their phenetic affinities, 
then it follows that as many taxa as possible 
should be treated in the analysis. I have used all 
the data available to me, and only H. palliatus 
galapagensis, H. p. pitanay and H. p. durnfordi 
are not represented. As some of the taxa used in 
this study are actually represented by individuals, 
I have followed Sneath & Sokal (1973) hereafter 
in employing the terminology operational taxo¬ 
nomic units (OTUs). 

The OTUs used in this study, the number of 
species in each and their geographic range in the 
breeding season are given in table I. The nomen¬ 
clature used is that of current checklists and con¬ 
notes no taxonomic decisions on my part. For 
heuristic purposes, and to provide a yardstick by 
which distances among OTUs can be judged, I 
have included hybridizing forms as separate OTUs. 
For example, the melanistic H. u. unicolor and the 
pied H. u. reischeki hybridize so extensively in 
northern New Zealand that they might well be 
lumped together in a single taxon if the extent 
of gene exchange is the final arbiter. Conversely, 
although H. bachmani hybridizes with H. palliatus 
in lower California, the hybrid population (here 
referred to as H. p. frazeri) is relatively restricted 
in distribution and thus the integrity of the gene 
pools of the parental species is essentially main¬ 
tained. 

Selection of characters 

Although numerical taxonomists such as Sneath & 
Sokal (1973) advocate the use of large suites of 
characters (commonly 60 or more) so that simi¬ 
larity coefficients will have reasonable confidence 
levels, in this study I have used only the seven 
external morphological characters traditionally 
used by avian systematists in assessing the affini¬ 
ties of closely related taxa. While I acknowledge 
that the use of a small number of characters may 
decrease the reliability of the statistical results, my 
intention is not to construct a numerical taxonomy 
as such, but rather to obtain a more complete eval¬ 
uation of the morphometric information contained 
in a traditional character suite. The characters used 
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Table I 


OTUs used in this study, their sample sizes and geographic ranges. 


OTU 

no. 

OTU designation 

Abbreviation 

9 9 

n 

3 $ 

Geographic range in the breeding season x ) 

i 

Haematopus ostralegus malacophaga 

H.o.mala. 

18 

17 

Iceland and Faeroe Islands 

2 

H. 

o. occidentalis 

H.o.occi. 

22 

21 

British Isles 

3 

H. 

o. ostralegus 

H.o.ostr. 

40 

37 

Europe east to Arkhangelsk and south to Spain 
and Baltic Sea; coasts of Asia Minor, Mace¬ 
donia, the Black and Caspian Seas 

4 

H. 

o. longipes 

H.o long. 

10 

1 

Inland southern and eastern Russia and western 
Siberia 

5 

H. 

o. finschi 

H.o.fins. 

18 

20 

South Island of New Zealand 

6 

H. 

o. osculant 

H.o.oscu. 

5 

11 

Coasts of northeastern Russia from eastern 
Siberia to northern China 

7 

H. 

o. longirostris 

H.o.rost. 

41 

45 

New Guinea south through Australia to Tas¬ 
mania 

8 

H. 

chathamensis 

H.chat. 

14 

14 

Chatham Islands 

9 

H. 

fuliginosus fuliginosus 

H.f.fuli. 

21 

28 

Coasts of Australia except ophthalmicus range 

10 

H. 

f. ophthalmicus 

H.f.opht. 

3 

4 

Northern Australia coast from Cape York 
Peninsula west to Kimberleys 

11 

H. 

unicolor unicolor 

H.u.unic. 

20 

21 

New Zealand 

12 

H. 

u. reischeki 

H.u.reis. 

17 

17 

North Island of New Zealand 

13 

H. 

moquini moquini 

H.m.moqu. 

8 

8 

Coasts of South Africa from Natal to South¬ 
west Africa 

14 

H. 

m. meadewaldoi 

H.m.mead. 

1 

1 

Eastern Canary Islands (probably extinct) 

15 

H. 

palliatus palliatus 

H.p.pall. 

4 

11 

Atlantic and Pacific coasts of North and 
Central South America 

16 

H. 

p. frazeri 

H.p.fraz. 

11 

18 

Lower California south to Jalisco 

17 

H. 

ater 

H.ater 

14 

29 

Coasts of South America from about 10°S on 
the west to about 43°S on the east; Falkland 
Islands 

18 

H. 

hachmani 

H.bach. 

17 

15 

Aleutian Islands through Alaska west to Lower 
California 

19 

H. 

leucopodus 

H leuc. 

Totals 

1 

285 

3 

321 

Southern South Africa from Chiloe Island 
—round to Chubut; Falkland Islands 


x ) These ranges are mapped more exactly in Baker (1973). 


and their methods of measurement are outlined 
below: (1) Bill length — measured along the 
dorsal surface of the oilmen from the tip to the 
junction of the feathers. (2) Bill depth — mea¬ 
sured from the junction of the feathers with the 
oilmen on the upper mandible to the exposed 
portions of the mandibular rami on the lower 
mandible. (3) Bill width — measured at the 
widest point of the bill beneath the dorsal feather 
margin on the upper mandible. (4) Tarsus length 
— measured diagonally from the posterior aspect 
of the joint between the tibiotarsus and the tarso- 
metatarsus, to the anterior aspect of the joint 
between the tarsometatarsus and the middle toe. 
(5) Middle toe length — measured dorsally from 
the joint with the tarsometatarsus to the base of 


the claw. (6) Tail length — measured from the 
point of emergence of the two central rectrices to 
the tip of the longest rectrix. (7) Wing length — 
measured as the chord of the flattened wing, from 
the carpal flexure to the tip of the outermost pri¬ 
mary feather. The first five measurements were 
taken with dial calipers to the nearest 0.1 mm, and 
the remaining two were measured with a steel 
rule accurate to 1 mm. 

Statistical methods 

As oystercatchers are markedly sexually dimorphic 
in size (Baker, 1975), character variation in males 
and females was analyzed separately beyond. The 
raw data matrix for each sex was too large to 
handle practically at the level of individuals, so 
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means were calculated for each OTU (except for 
those where n = 1) and the dimensions of the 
data matrix were correspondingly reduced to 
19 X 7 (number of OTUs X number of char¬ 
acters) . 

Prior to the calculation of similarity coefficients, 
the characters were standardized across OTUs so 
that they had a mean of zero and a standard 
deviation of one, and thus each character theoreti¬ 
cally can contribute equal information in the en¬ 
suing analyses. Average taxonomic distances (d) 
and product-moment correlation coefficients (r) 
were then computed among OTUs, and separate 
cluster analyses were performed using the un¬ 
weighted pair-group method on arithmetic aver¬ 
ages (UPGMA — Sneath & Sokal, 1973). 
UPGMA clustering was initially used because it 
usually gives the least distortion when represented 
by a two-dimensional phenogram (Sokal & Rohlf, 
1962; Farris, 1969; Rohlf, 1970). Cophenetic 
correlations between the phenograms and the orig¬ 
inal distance and correlation matrices were cal¬ 
culated to determine how accurately the pheno¬ 
grams summarized the similarity coefficients. 

The results of hierarchical clustering were tested 
against a non-hierarchical ordination technique 
combining both principal component analysis 
(PCA) and nonmetric multidimensional scaling 
analysis (MDSCALE). This dual approach utilizes 
the results of PCA as an initial input configuration 
for the MDSCALE routine. 

The value of principal component transforma¬ 
tion of morphometric data derives from the assess¬ 
ment of simultaneous variation in all characters 
over all OTUs, and the consequent reduction of 
a large data matrix to a much smaller set of 
orthogonal (and thus uncorrelated) principal com¬ 
ponents. The MDSCALE technique, on the other 
hand, represents t OTUs geometrically by t points 
in a reduced space of k dimensions (where k < t; 
k typically = 2 or 3) such that the interpoint 
distances have a monotone relationship with the 
distances among the OTUs in the original dis¬ 
similarity matrix. MDSCALE differs from PCA 
in requiring only that the monotone relationship is 
preserved in the reduced space, rather than ex¬ 
plaining a maximum percentage of the variation. 
Results of the two techniques are in practice quite 


similar, but MDSCALE more faithfully preserves 
the distances between the OTUs in the reduced 
space (Rohlf, 1970). The goodness of fit of the 
MDSCALE solution to the distance matrix was 
assessed by computing stress values (see Kruskal, 
1964) and the matrix correlation r M * (Rohlf, 
1972). 

In this study the OTUs were projected into a 
three-dimensional space, and a minimum spanning 
tree (MST) computed from the original distance 
matrix was superimposed to provide additional 
insights on the phenetic relationships of the OTUs. 
Distinctive clusters of OTUs in the three-dimen¬ 
sional MDSCALE ordinations were isolated using 
subsets analysis. This technique finds all subsets 
of OTUs satisfying the condition that the greatest 
dissimilarity between OTUs in the subset is less 
than the least dissimilarity between any OTU in 
the subset and any OTU not in the subset. Only 
the complete set of OTUs and single OTUs are 
not permitted as subsets. 

RESULTS 
Cluster analysis 

The morphometric relationships of the 19 OTUs, 
as elucidated by UPGMA cluster analyses of the 
r and d matrices, are shown for adult females and 
males in fig. 1. Results for both sexes are quite 
disparate for the two coefficients of similarity. 

For females, three clusters are produced from 
correlation coefficients and two from average 
taxonomic distances. Moreover, the composition 
of clusters changes somewhat for the two coeffi¬ 
cients. While cluster analysis of the r matrix (fig. 
la) faithfully preserves the known close phenetic 
relationships between the hybridizing forms of 
northern New Zealand ( H. u. unicolor and H. u. 
reischeki) and lower California (H. p. palliatus 
and H. p. jrazeri), it has the disruptive effect of 
separating the OTUs representing the western 
European subspecies of H. ostralegus into different 
clusters. More specifically, it breaks off H. o. 
occidentalis and H. o. malacophaga (British, Ice¬ 
landic and Faeroe birds) from the cluster including 
the nominate H. o. ostralegus (the remaining 
western European birds). The central European 
and Eurasian subspecies (H. o. longipes and H. o. 
osculans respectively) group with H. o. ostralegus. 
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Fig. 1. Phenograms from UPGMA cluster analysis of similarity matrices: (a) correlation phenogram for females, (b) distance 
phenogram for females, (c) correlation phenogram for males, (d) distance phenogram for males. OTUs are identified as in 
table I and cophenetic correlations are denoted by r c . 


and with H. o. jinschi of New Zealand. The 
Australian Pied Oystercatcher ( H. o. longirostris ) 
clusters with western European H. ostralegus and 
the melanistic H. m. moquini of South Africa. 
Both Australian subspecies of the Sooty Oyster- 
catcher ( H. fuliginosus) group together. The Chat¬ 
ham Islands Oystercatcher ( H. chathamensis) 
appears as an aberrant form distantly related to the 
Canary Islands H. m. meadewaldoi. 

Results of the cluster analysis of the matrix of 
average taxonomic distances (fig. lb) are intu¬ 
itively more appealing (though note the comments 
of Rohlf & Sokal (1965) to the contrary regarding 
size-ordered classification from d matrices), not 


only because they are largely consistent with taxon¬ 
omies adopted by recent checklists (e.g. Kinsky, 
1970; Condon, 1975), but also because they are 
explicable in terms of zoogeographic and evo¬ 
lutionary processes. All the subspecies of H. ostra¬ 
legus except H. o. longirostris fall within a major 
cluster, and the New Zealand H. o. jinschi shows 
closest phenetic affinities with the Eurasian H. o. 
longipes and H. o. osculans. H. chathamensis of 
the Chatham Islands clusters with western Euro¬ 
pean subspecies of H. ostralegus and the New 
World H. leucopodus and H. bachmani. Both 
Australian subspecies of H. fuliginosus group 
together and link with the South American Black 
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Oystercatcher (H. ater) and the Australian Pied 
Oystercatcher ( H . o. longirostris). The hybridizing 
forms of both H. unicolor and H. palliatus form 
compact groups and together form a broader 
cluster, with the South African H. m. moquini 
and the Canary Islands H. m. meadewaldoi as 
more distant outliers. 

Neither of the phenograms are particularly good 
in summarizing the information in the similarity 
matrices. Cophenetic correlations are only 0.735 
for the correlation coefficient phenogram and 
0.743 for the distance phenogram. Cluster analysis 
of OTUs of males, as with females, produces three 
clusters with correlation coefficients and two with 
average taxonomic distances. The membership of 
clusters differs between the sexes, but the overall 
patterns of groupings of OTUs are broadly similar. 
Cluster analysis of the r matrix again partitions 
the subspecies of H. ostralegus into different 
clusters, but more disturbingly H. chathamensis 
is inserted between the hybridizing forms of H. 
unicolor (fig. lc). The integrity of the H. palliatus 
complex is maintained, however, and the two sub¬ 
species of H. fuliginosus group together. Western 
European subspecies of H. ostralegus appear to be 
phenetically closest to the South African H. m. 
moquini and the Australian H. o. longirostris. 

The results of cluster analysis on the d matrix 
of males (fig. Id) are again intuitively much more 
satisfactory than those from the r matrix. All the 
Eurasian subspecies of H. ostralegus except H. o. 
osculans appear in the same cluster, along with the 
ostralegus -like H. leucopodus of South America 
and H. chathamensis of the Chatham Islands. H. o. 
osculans is a distant outlier in a cluster including 
the hybridizing forms of H. palliatus from lower 
California and the Canary Island Oystercatcher 
(H. m. meadewaldoi). The OTUs representing 
H. unicolor and H. fuliginosus are clustered to¬ 
gether with H. ater of South America. 

Both phenograms of males provide relatively 
poor fits with the similarity values in the r and d 
matrices. The cophenetic correlations are only 
0.688 for the correlation coefficients and 0.695 
for the average taxonomic distances. This dis¬ 
appointing result reflects the weakly hierarchical 
nature of the similarity matrices, which in turn 
derives from the relatively constrained and homo¬ 


geneous range of phenetic variation found in the 
Haematopodidae. 

Principal components analysis 

Because the PC A projections of the OTUs are 
quite similar to those of MDSCALE in this study, 
they will not be presented here. However, it is 
useful to inspect the pattern of loadings of the 
characters on the first three principal components, 
as they aid considerably in interpreting the axes 
of the MDSCALE ordination. For females, 91.7% 
of the character variance is accounted for on the 
first three principal components; 59.8% on com¬ 
ponent I, 20.2% on II and 11.7% on III. For 
males, 86.9% of the character variance is ac¬ 
counted for by the first three principal com¬ 
ponents; 54.4% on component I, 21.2% on II and 
11.3% on III. 

The loading patterns are similar for the sexes 
(table II). Component I has high positive loadings 
for all characters except bill length, and therefore 
serves to differentiate the OTUs largely on the 
basis of size. Component II is effective in dis¬ 
tinguishing OTUs chiefly on an inverse relation¬ 
ship between bill length (and tarsus length to a 
lesser degree) and tail length. Middle toe length 
and bill width load positively, and wing and tail 
lengths load negatively on component III of fe¬ 
males, and the reverse is apparent for males. Hence 
component III seems to involve similar contrasts 
in these characters across OTUs. 

Table II 


Loading of characters on the first three principal components 
based on a matrix of correlations among morphological 
characters. 


Character 

Sex 

Com¬ 
ponent I 

Com¬ 
ponent II 

Com¬ 
ponent III 

Bill length 

Male 

0.163 

—0.880 

0.405 


Female 

0.236 

—0.919 

—0.261 

Bill depth 

Male 

0.849 

—0.248 

—0.031 


Female 

0.894 

—0.142 

—0.103 

Bill width 

Male 

0.837 

—0.220 

—0.381 


Female 

0.901 

—0.088 

0.320 

Tarsus length 

Male 

0.804 

0.304 

0.132 


Female 

0.808 

—0.411 

0.056 

Middle toe length 

Male 

0.797 

0.140 

—0.433 


Female 

0.763 

0.184 

0.575 

Wing length 

Male 

0.821 

0.350 

0.296 


Female 

0.863 

0.277 

—0.376 

Tail length 

Male 

0.634 

0.605 

0.433 


Female 

0.733 

0.512 

0.402 
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Nonmetric multidimensional scaling and subsets 
analysis 

The three-dimensional MDSCALE ordinations of 
females and males differ somewhat in the place¬ 
ment of particular OTUs, but the general pattern 
of phenetic relationships expressed therein is 
similar (fig. 2). In both ordinations the Eurasian 
subspecies of H. ostralegus group together and 
have clear phenetic affinities with H. o. finschi 
of New Zealand, H. leucopodus of South America 
and H. bachmani of western North America. All 
the remaining Australasian taxa group together 
within another major grouping containing all the 
New World OTUs. The African H. moquini 
is phenetically intermediate, appearing nearer the 
New World group of females (fig. 2a) and the 
Old World group of males (fig. 2b). 

Distinctive clusters of OTUs within the two 
major groups are indicated by the subsets proce¬ 
dure. For females, the western European sub¬ 
species H. o. occidentalis and H. o. malacophaga 
are very similar phenetically, as are the Eurasian 
subspecies H. o. longipes and H. o. osculans. As 
might be expected, for both sexes the New World 
and New Zealand hybridizing forms fall into 
specific subsets, as do the subspecies of the Austra¬ 
lian H. fuliginosus. A slightly different group of 
subsets emerges for Old World males; H. o. occi¬ 
dentalis has closest phenetic affinity with H. o. 
ostralegus, and H. o. finschi of New Zealand 
appears in the same subset as H. leucopodus of 
South America. 

The low stress values (s) and the high matrix 
correlations (rm*) indicate that the ordinations 
provide an excellent fit to the phenetic relation¬ 
ships represented in the matrix of average taxo¬ 
nomic distances. Given the small sample sizes on 
which some OTUs are based, and the small number 
of characters used, the degree of concordance of 
the ordination patterns of males and females at¬ 
tests to the robustness and utility of the joint PCA/ 
MDSCALE approach. 

DISCUSSION 

Taxonomists have long been impressed by the 
distinctive but uniform morphology of oyster- 


catchers, and have accordingly classified them in 
a single monogeneric family. Consensus has ended 
at the generic level, however, because of the 
seemingly intractable nature of morphometric 
variation of allopatric forms, the lack of a global 
view in assessing affinities, and the confounding 
influence of melanism which has probably arisen 
independently in several geographically wide¬ 
spread taxa (e.g. H. unicolor in New Zealand, 
H. fuliginosus in Australia, H. moquini in Africa, 
H. ater in South America and H. bachmani in 
North America). 

Taxonomic problems within the Haematopodi- 
dae are thus similar to those encountered within 
homogeneous assemblages such as species-com¬ 
plexes, and their successful resolution often de¬ 
pends on a careful synthesis of relevant data from 
morphology, ecology, ethology, physiology etc. 
(Mayr, 1969). This approach has been fruitful 
for the New Zealand taxa (Baker, 1974), but 
reliable data from such diverse disciplines as these 
are as yet unavailable even for the Australasian 
taxa. 

Yet, there is still much valuable taxonomic in¬ 
formation to be gleaned from morphological data 
sets. The use of non-hierarchical ordination 
methods seems particularly appropriate for analyz¬ 
ing phenetic relationships of homogenous taxa 
as the methods do not force them into a rigid 
and possibly artificial hierarchy. In this study 
multivariate analysis of an admittedly small set of 
characters long used by avian systematists has given 
new insights into the nature of morphological 
variation in oystercatchers. Of particular impor¬ 
tance in this respect is the finding that two principal 
groups of OTUs can be defined, one an Old 
World group with a solid nucleus of Eurasian 
H. ostralegus and the other a New World group, 
which also includes all the Australasian OTUs 
except H. o. finschi of New Zealand. This latter 
taxon seems better considered as a Palaearctic 
element in the New Zealand fauna (following a 
secondary invasion) rather than a Holarctic ele¬ 
ment as suggested by Fleming (1962). Conversely, 
the Australian H. o. longirostris appears to be 
phenetically similar to the New World forms, and 
should probably be given full species rank. How¬ 
ever, zoogeographic and taxonomic speculation 
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Fig. 2. Three-dimensional ordinations using the PCA/MDSCALE method: (a) females, (b) males. Dimension III is 
represented by the height of the projections. OTUs are identified as in table I, stress is denoted by s and the matrix 
correlation by r^*. The MST is indicated by the dashed lines, and joins phenetically similar OTUs together. Distinctive 
clusters of OTUs defined by the subsets procedure are enclosed by solid lines. 
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based on the present data set is obviously prema¬ 
ture, and must await the results of a study of a 
much larger set of skeletal and biochemical char¬ 
acters, currently in progress. 
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